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We employ three dimensional x -ray coherent diffraction imaging to map the lattice strain distribu- 
tion, and to probe the elastic properties of a single crystalline Ni (001) nanowire grown vertically 
on an amorphous Si02 || Si substrate. The reconstructed density maps show that with increasing 
wire width, the equilibrium compressive stress in the core region decreases sharply while the surface 
tensile strain increases, and gradually trends to a nonzero constant. We use the retrieved projection 
of lattice distortion to predict the Young's Modulus of the wire based on the elasticity theory. 
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Due to their high surface-to-volume ratio, transition 
metal nanostructures such as nanowires (NWs) could 
potentially be used in a broad range of applications in 
catalysis .sensors, batteries, fuel cells, and magnetic de- 
vices. [l|-[4| Most of these devices require NWs with well- 
defined size, shape, and spatial ordering. In addition, 
economical routes to mass production of NWs are desir- 
able for practical applications. However, structural mor- 
phology of the NWs is strongly affected by the complex 
interplay between several growth parameters, such as 
temperature of the substrate, vertical and lateral mate- 
rial transfer, growth rate, in-plane mobility of ad- atoms, 
etc. In order to understand resulting morphology of NWs 
it is therefore vital to develop characterization tools that 
can nondestructively probe the three dimensional (3D) 
structural and mechanical properties of these NWs with 
nanoscale resolution. 

In this letter we study the structure, shape, 3D lat- 
tice distortion, and the electron density distribution in 
single crystal Ni NWs with the aid of synchrotron-based 
coherent x-ray diffraction (CXD) imaging Q. The 
NWs were grown using thermal chemical vapor deposi- 
tion (CVD) on a Si substrate that has been oxidized re- 
sulting in a 500-nm Si02 amorphous coating layer. CVD 
at 650 °C yields densely populated coverage of vertically- 
oriented single-crystal NWs depicted in Fig. []Ib) with 
well-defined orthogonal and smooth facets. [3, Q The 
NWs are grown in the [001] crystallographic orientation 
with widths ranging from 50 nm to about 300 nm and 
lengths of up to 5 pm and average coverage density of 
between 0.1 and 0.3 NW/pm 2 . 

CXD in the Bragg geometry has been shown to be 
a powerful characterization technique for imaging of lo- 
cal nanoscale lattice distortions, commonly described by 
deformation fields, u(r) within a small crystal. ifl 
CXD allows us to image not only the overall shape of the 
nanostructure in 3D, but also the projection of the crystal 
lattice displacement field on to the Q vector of the mea- 
sured Bragg spot. The lattice distortion represents the 
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imaginary component of object density and manifests it- 
self as local asymmetry of the measured Bragg diffraction 
pattern. Upon an inversion of coherent x-ray diffraction 
pattern, the resulting real-space image of the object will 
be complex valued p(r) = pG hkl (r) exp [-iG hk i • u(r)], 
where the amplitude pG hkl represents the density of the 
object and the phase, Ghki • u (r) is the projection of the 
atomic displacement vector field u (r) relative to atomic 
positions in an ideal perfectly periodic lattice onto the 
reciprocal lattice vector Ghki for the measured Bragg re- 
flection (hkl). Since only the intensities, and not the 
phases, are measured experimentally, direct inversion of 
the CXD pattern to obtain a real-space distribution of 
the object density is impossible due to phase problem. 
However, the phase problem [9] can be solved via spe- 
cial phasing procedures, provided that the scattering ob- 
ject is finite and that diffraction pattern is sufficiently 
oversampled. 0, 0, 

Figure QJa) shows the experimental setup used for 
CXD at Beamline 34-ID-C at the Advanced Photon 
Source (Argonne National Laboratory). A monochroma- 
tor was used to select E = 8.919 keV x-rays with 1 eV 
bandwidth, resulting in a longitudinal coherence length 
of about 0.1pm. Slits were placed into the beam up- 
stream of the Kirkpatrick-Baez focusing mirrors to select 
a 20pm x 50pm wide transversely coherent segment of 
the beam that was focused down to about 1.5pm to illu- 
minate a single NW. Coherent diffraction patterns were 
recorded for the rocking curves of the (004) and (111) 
Bragg reflections by rotating the sample through the 
Bragg condition in increments of about 0.005 degrees. 
A CCD detector with 20pm pixel size was used to col- 
lect the 2D diffraction slice for each rotation angle. Full 
3D diffraction patterns were then constructed by tomo- 
graphically stacking these 2D frames together. 

The lack of inversion symmetry for CXD patterns such 
as the one shown in Fig Ufa) and Fig J2] are a clear indica- 
tion that the imaged nanostructure is strained, resulting 
in non-zero imaginary component of effective density. [14| 
We inverted our CXD patterns using a combination of 
Fienup's Hybrid Input-Output [i| and Error Reduction 
algorithms. Electron diffraction patterns (see Fig. \Hb)) 



FIG. 1. (Color online) Schematic of the experimental set-up and sample, (a) Monochromatic X-ray beam (X) of wave vector 
kin impinges on to Kirkpatrick-Baez mirrors (not shown on sketch), which creates a localized illumination on the sample (S). 
By rotating the sample through the Bragg condition in increments of about 0.005 degrees, CXD patterns in the vicinity of the 
(111) reciprocal lattice point are recorded with a two-dimensional pixelated detector D. A typical diffraction pattern (out of the 
hundreds of patterns collected) is shown here, (b) Selected area electron diffraction pattern from a vertical Ni NW indicates 
that it is a single crystal with a 001 oriented growth axis, (c) and (d) Scanning electron micrographs of the sample reveal that 
the NWs are well oriented with smooth facets, (e) Basal part of the NW exhibits pronounced tapering features and orthogonal 
faceting. 



and SEM (Fig. Hfc-d)) indicate that the NWs are defect- 
free, vertically oriented (grown along 001) single crys- 
tals with smooth facets. A priori knowledge of the NWs 
shapes from SEM and the thickness oscillations along the 
vertical direction clearly seen from the diffraction pat- 
terns in Fig.QJa) and Fig. [2] were utilized to determine a 
good starting guess of the size of the region (support) at 
which the reconstructed object is allowed to exist. 

From a priori NWs growth information, @, [1| we ob- 
served that the NWs had a small amount of impurities 
embedded at the surface. Taking this into account, we 
introduce an inner core region of the NW to be of uni- 
form density and having relatively few impurities while 
the outer shell region has a higher density of impuri- 
ties. The x-ray diffraction patterns from this core-shell 
model are composed of an interference sum of individ- 
ual wave field contributions scattered from the two re- 
gions: core and shell of the NW. We first define an un- 
perturbed state of the crystal by assuming that for each 
region there is idealized homogeneous composition and 
strain, spatially limited by the shapes of the individual 
crystals. We then introduce the perturbation by taking 
into account the effects of non-uniform strain relaxation 
and possible compositional fluctuations within the pre- 
viously homogeneous region of the crystals. The contri- 
bution of each wave field is the result of the convolution 
of two Fourier transforms: (1) the individual reciprocal 
shape function which leads to intensity decay of at least 
q~ 2 with respect to the unperturbed Bragg peak of that 
material and (2) the Fourier transform of p (r) contain- 
ing the spatial fluctuations in amplitude (compositional 
fluctuation) and phase (displacement field). 



The spread of the latter Fourier transform in reciprocal 
space (see Fig. [2]) is dependent on the maximum values of 
the spatial gradients of the displacement field in the ide- 
alized layers and is related to the strain-induced broad- 
ening of the diffracted signal AQ P (along a given p direc- 
tion) near the Bragg spot. A more detailed description of 
this approach can be found in our earlier work. [l3|, [l5| 
We transform this band-limited frequency domain con- 
straint into a direct space constraint of limited lattice 
displacement gradient. By utilizing the uniformity of the 
electron density (and its threshold) constraints [15] along 
with the standard support constraints, we reconstruct the 
3D distribution of NW electron density (shape) and cor- 
responding phases (related to the projection of displace- 
ment field) from the (004) and (111) Bragg peaks, shown 
in Figs. QJi). We estimate the magnitude of the maxi- 
mum values of displacement derivative along the vertical 
direction (from 004 Bragg peak) and along the (111) di- 
rection (from 001 Bragg peak). These values were used 
as a limit for the maximum allowable phase difference 
in real space for neighboring points in the axial and the 
(111) directions. 

Figures [3fa) and[3fb) show the refined support (outer 
shape) and the reconstructed phases of the NW while the 
reconstructed amplitudes (internal distribution of elec- 
tron density) are displayed in Fig. [3fc). The solution to 
our phasing problem puts the NW width at about 220 nm 
and length of 2 /im, consistent with SEM measurements. 
A 3D reconstruction of the distribution of phases from 
the (004) and (111) Bragg peaks are shown in Fig. SJa). 
Unwrapping these phase distributions gives us the mag- 
nitude of the axial (001) and radial (110) components 
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FIG. 2. (Color online) Measured reciprocal space map (log 
10 of intensities) from Ni NW near the 004 Bragg reflection. 



of the displacement field Fig. Hfb). The color-bar shows 
that the maximum average displacement of 0.15 nm (0.4 
of the Ni lattice parameter) is attained at the top of the 
NW. The fact that the bottom of this NW is contracted 
may be attributed to the formation of free surfaces from 
the NW tapering. The height of the tapered section will 
depend on the specific sample, but the angle between the 
main body and the taperingportion could be formed by 
the (110) and (100) planes. 0, Is| To test this assumption, 
we utilized the Bragg peak broadening in the vicinity of 
(111) and (004) to obtain a measure of the maximum 
displacement gradient along those crystallographic direc- 
tions. [ISElI 

In order to confirm the reproducibility and uniqueness 
of the obtained solutions, we performed a series of phase 
retrieval procedures with different input random set of 
phases for the measured scattered radiation. We defined 
a cost function in the form of an error metric as E\ = 

Eti(l*f m l " V / ^) 2 /Ef =1 1? XP , where \Ff m \ is the 
magnitude of the simulated amplitude and lf xp is the 
experimental intensity of point i in the reciprocal space 
map. 

The phasing process, for a given initialization of the 
algorithm, was determined to be complete when the er- 
ror metric is reduced to 10 -8 . The maximum number of 
iterations allowed was typically about 4,000. The differ- 
ences between many equivalent solutions obtained from 
random initial sets phases provide a measure of the reso- 
lution of the resulting NW reconstruction in the form of 
a Phase Retrieval Transfer Function. [12j The real-space 
resolution attained in this work is about 20 x 20 nm 2 . 

From Figs. HJc) and HJd) we notice the presence of a 
tensile strain at the near-surface regions and at the tip of 
the NW, contrasted with significant compressive strain 
in the core. We observe that the tensile strain decays 
non-uniformly along the z-axis while changing smoothly 
in the core region to a compressive strain at the bot- 
tom of the NW. At the top-shell region of the NW (see 
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FIG. 3. (Color online) 3D reconstructions of (a) support, (b ) 
phases and (c) amplitudes of the NW complex valued density 
function from the (004) Bragg peak. 

Fig. ID^d)), the strain is purely tensile but nonzero. We 
also notice (see Fig.HJc)) that the core region of the NW 
has a resultant uniaxial compressive lateral displacement 
field ranging from zero to about -0.15 nm while the max- 
imum tensile displacement magnitude is 0.15 nm. The 
core region of the NW shows a clear homogenous density 
distribution while at the shell it varies slightly, as shown 
in Fig. [3fc). This can be attributed to the presence of 
surface impurities (likely Si) during NW growth. These 
impurities are associated with a slight change in the aver- 
age electron density. This causes a change of the optical 
path length of the x rays and thereby a phase variation. 
However, this phase information alone is insufficient to 
draw definitive conclusions on the quantitative impurity 
contents. 

From the retrieved displacements, utilizing the elastic- 
ity theory [16| and experimentally verified values for the 
stiffness constants, we obtain quantitative change in the 
stress from —3.6 GPa at the center of the NW to about 
13.4 GPa in the outermost part. However, the stress re- 
mains approximately constant within a core region of the 
NW where the density is highly uniform. This nonzero 
stress is due to residual strain in the NW.. This implies 
that the surface endures the tensile strain while the inte- 
rior is compressive, even if this NW is in equilibrium, thus 
providing a possible explanation for observed tapering. 

Most metals have the possibility of exhibiting nonlin- 
ear and yet elastic behavior i.e., the stress a is essentially 
nonlinear to strain within the elastic limit. [l6| Under 
such conditions, the stress-strain relation can be approx- 
imated by a Taylor series expansion a = ao + aiS + a2e 2 + 
(e 3 ) ~ a±e + ci2£ 2 . The term <jo is negligible if the 
crystal is unstrained. For extremely small strain, i.e., 
if the length-scale over which the deformation varies is 
much larger than the discrete length of the matter [l6| 
(in this case much larger than the lattice parameter of 
Ni), the stress functions obey a quadratic strain behav- 
ior. By adopting this quadratic type variation, we obtain 
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the quantitative value of the Young's modulus of Ni NW 
given by a\ to be 100 GPa while its dependence due to 
strain given by a2 is 480 GPa. 



In conclusion, we have used CXD to investigate the 
shape and internal structure of CVD-grown single crys- 
tal Ni NW grown on an amorphous SiC>2 substrate. We 
demonstrate that the basal NW taper is accompanied by 
substantial residual strain within the NW, as evidenced 
by a transition from compressive strain in the core of 
the NW to tensile strain in the near-surface shell. The 
quantitative values of the elastic properties such as the 
Young's Modulus and its dependence on strain have been 
calculated from the imaged strain field using the elas- 
ticity theory. Characterization of the Ni NWs confirms 
their exceptional crystalline quality and indicates their 
suitability for numerous wide-ranging applications, while 
strain engineering opens a gateway to novel NW-based 
device functionalities under external stimuli such as pres- 
sure, electric or magnetic fields. 



FIG. 4. (Color online) 3D reconstructions of (a) the phases, 
and (b) the un- wrapped phases from (004) and (111) Bragg 
peaks for a single crystal Ni NW converted to projection of 
the vertical and axial components of the displacement field, 
(c) The variation of vertical component of displacement field 
along the axial (001) direction within the core (blue) and the 
shell (red) of the NW. (d) The variation of the vertical com- 
ponent of displacement field along the radial (110) direction 
of NW width at the top (purple) , middle (black) and bottom 
(green) parts of the NW, at select regions indicated in b). 
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